Design of a far-infrared spectrometer for atmospheric 
thermal emission measurements 


David G. Johnson 

NASA Langley Research Center, Hampton, Virginia, USA 


ABSTRACT 

Global measurements of far infrared emission from the upper troposphere are required to test models of cloud radiative forcing, 
water vapor continuum emission, and cooling rates. Spectra with adequate resolution can also be used for retrieving atmo- 
spheric temperature and humidity profiles, and yet there are few spectrally resolved measurements of outgoing longwave flux at 
wavelengths longer than 16 gm. It has been difficult to make measurements in the far infrared due to the need for liquid-helium 
cooled detectors and large optics to achieve adequate sensitivity and bandwidth. We review design considerations for infrared 
Fourier transform spectrometers, including the dependence of system performance on basic system parameters, and discuss the 
prospects for achieving useful sensitivity from a satellite platform with a lightweight spectrometer using uncooled detectors. 
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1. INTRODUCTION 

Fourier transform spectrometers (FTS) operating in the far infrared (FIR; defined here as wavelengths longer than 16 gm) have 
generally used cooled detectors, either photoconductors operating at 4.2 K or bolometers operating at even colder temperatures. 
The low temperatures are required both to achieve useful sensitivity and, in the case of a scanning FTS, to provide the neces- 
sary electronic bandwidth. Recent improvements in the sensitivity and bandwidth of uncooled detectors have motivated us to 
reconsider the possibility of demonstrating useful sensitivity with an FTS operating in the FIR with uncooled detectors. 

In the following sections we derive a general expression for the sensitivity of a Fourier transform spectrometer as a function 
of basic system parameters. We then estimate the sensitivity that can be achieved with the present generation of uncooled 
detectors, and the improvement in performance that could be realized with detectors operating at the thermal noise limit. 

2. INTERFEROGRAM 

We consider first a conventional scanning FTS, shown schematically in Figure 1 . For a monochromatic source the interferometer 
output is given by 


I(x) =Io + 0.5q(k)SrfCos(27tLv) , (1) 

where I(x) is the output power, x is the optical path difference (OPD) between the two arms of the interferometer, Io = S 2 + 
0.5r|(k)S</, S f / = Si — S 2 , T| (k) = 4 R(k)T(k), Sj and So are the power at each input, k is the wavenumber (1/wavelength), and 
R(k) and T (k) are the beamsplitter reflectance and transmittance, respectively. 

For a broadband source the interferogram is given by the integral of Equation (1) as shown below: 

/ +°° 

dki\(k)Sd(k)cos(2nkx ) , (2) 

/ +“ 

d*[S 2 (*)+0.5Ti(*)S (/ (*)]; (3) 
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where Sd(k ), the net input spectral power per unit wavenumber interval, is given by Si(k) — 5o(k), and Sj (k) and 5o(k) are the 
spectral power per unit wavenumber interval at inputs 1 and 2, respectively. 
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Figure 1. Generic Fourier transform spectrometer layout. 

We estimate the net input spectrum by taking the Fourier transform of Equation (2) after subtracting the constant term 7o: 

/ +oo f+°° 

dxexp(2nikx) / dk'0.5’r\(k!)S d (k!)cos(2nk!x) , (4) 

= 0.5r }(k)S d (k); (5) 

where S' (k) is the estimated net input spectrum. 

In practice the inputs Si and S 2 (or the output I(x)) are optically filtered to produce a band-limited interferogram, and we 
further reduce the noise by passing the detector output through an electronic filter matched to the optical bandpass and mirror 
scan velocity. We then estimate S' c/ (k) from the discrete Fourier transform (DFT) of the band-limited interferogram 1 : 

S' d (k n ) ~ AH „ ; (6) 

N/ 2-1 

H n = Y, exp(27tk„x j) (Ij — Iq) , (7) 

j=-N/2 

where H n is the DFT of (Ij —Io), Ij = /(xy), Xj = /A, k„ = n/NA , N is the total number of samples, and A is the difference in 
optical path between interferogram samples. Note that we have assumed that the interferogram is fully two-sided. 

We combine Eqs. (5) and (6) to solve for //„ as a function of the net input power per unit wavenumber: 

H n &0.5r\(k)S d (k„)/A. ( 8 ) 


To estimate the response of a particular interferometer to an input of known net spectral intensity we need to consider 
optical throughput and spectral response of the detector. To simplify the present discussion we assume that the spectrometer 
fully illuminates a detector of area A with a beam having the focal ratio /, and that the spectral response of the detector, filters, 
and other optics (except the beamsplitter) is represented by a single efficiency e(k). The effective (detected) net spectral power 
per unit wavenumber ( S d ) is given by: 


S d (k) = E(k)F d (k)nA/4f 2 , 


where F d is the net input spectral intensity. 



3. NOISE 


We now consider the effect of noise on the interferogram. We assume that noise is given by the function g(x), and that the 
observed interferogram is given by I(x) + g(x). We define gj = g(xj ), where Xj is defined above. The variance in the DFT of 
gj is related to the variance of gj by the discrete version of Parseval’s Theorem : 

Var(G„)=AWar(g;), (10) 

where G n is the DFT of gj, Var(y) signifies the variance of y, and N is defined above. The variance of gj is related to the specific 
detectivity and electronic bandwidth of the detector and preamplifier as follows: 

Var ( gj )=Av B /(D*) 2 , (11) 

where Vg is the preamplifier bandwidth and D is the specific detectivity of the detector. After substituting Equation (11) in 
Equation (10) and taking the square root we derive the result 

o = s/NAvb/D* , (12) 

where a = \/Var(G„). Note that both G„ and H n are complex (having both real and imaginary components), while g j and 
the original spectrum S d are real. We can reduce o by an additional factor of \/2 by properly phase correcting the spectrum 
(calculating 9 X{H„ exp [/(])(/:„)] }, where the phase function </>(k„) can be estimated in a number of ways 2 ), since half the variance 
in G„ corresponds to the imaginary component of .S',/ and can be removed. 

The signal-to-noise ratio (SNR) is given by H„/a. Substituting Eqs. (8) and (12) for H n and o, respectively, we derive the 
result 


SNR = 0.5D*j\(k)S d (k) /A^/NA vb ■ (13) 

Note that the unapodized spectral resolution of a transformed two-sided interferogram of length N is given by 1 /NA. We now 
assume that Vg is equal to the Nyquist frequency (meaning that we sample the interferogram at exactly 2 Vg) so that the time 
required to record one interferogram is given by /V/2vg. We define 8k = 1 /NA and x = N /2v/j. After combining Eqs. (9) and 
(13) and substituting for x and 8 k, we derive the result 

SNR = TK(k)D*^k)F d (k)8kV2%A/8f 2 . (14) 

4. SENSITIVITY 

The sensitivity (NEdP) is given by G//SNR, or 

NEdP = 8f 2 /mD*y\(k)8k\/2iA , (15) 

where NEdP is expressed in units of power per unit area per unit solid angle per unit wavenumber interval. The minimum 
detectable change in temperature is given by NEdP/(dF /dT), where dF /dT is the derivative with respect to temperature of 
the spectral intensity of a blackbody source. While NEdP depends only on the characteristics of the spectrometer and detector, 
NEdT also depends on the assumed source temperature. 

We now consider the steps that can be taken to maximize the sensitivity of an FTS operating in the FIR with uncooled 
detectors. To be more specific, we will consider an imaging FTS in low earth orbit. The instrument is to provide spectra 
covering the wavelength range 10-100 /rm with 0.6 cm -1 resolution (apodized), daily global coverage, and a 10 x 10 km 
footprint. The detectivity of uncooled far-infrared detectors depends on detector temperature and is independent of source 
intensity, 3 so that we may assume that each of the terms in Equation (15) is independent. 

The need for daily global coverage and a small footprint greatly limits the integration time x. Observing the entire globe with 
a 100 km 2 footprint once a day means that the spectrometer must provide five million spectra per day, corresponding to 17 ms 
per spectra if they are obtained one at a time. To obtain the desired resolution of 0.6 cm -1 (corresponding to 8 k = 0.3 cm -1 ) 
in 17 ms we must scan the interferometer mirror at 98 cm/s, resulting in an electronic frequency of 98 kHz at a wavelength 
of 10 gm. Available uncooled IR detectors simply do not provide adequate sensitivity at such high modulation frequencies. 
However, by going to a modest 10x10 array of detectors we increase the integration time and reduce the modulation frequency 
by a factor of 100 to a more reasonable 980 Hz. The penalty is that the FTS throughput must also be increased by a factor 




Figure 2. Comparison of conventional and bilayer pellicle beamsplitter. The conventional beamsplitter consists of an uncoated 
12 pm thick Mylar film, while the bilayer beamsplitter consists of 3.5 pm of polypropylene coated with 1.15 pm of germanium. 


of 100, although this is not an insurmountable problem. In the discussion that follows we assume that the FTS illuminates a 
10 x 10 array, giving us 1.7 s per field, and that of that 1.7 s we have 1.4 s available to record the interferogram. This gives 
X = 1 .4 s and Vg = 1200 Hz for a minimum wavelength of 10 pm. 

Clearly the detector needs to be illuminated with the fastest possible optical system since NEdP is proportional to f 2 . The 
best that can be done with a reflective off-axis system is about / = 2. While an on-axis / = 1 refractive system could be built, 
it is unlikely that such a system could be color-corrected over the desired optical band. However, one can do much better 
with non-imaging reflective systems. A parabolic feed horn 4 provides the maximum possible flux concentration, achieving 
/ = 0.5. Such feed horns are commonly used with liquid-helium cooled FIR detectors. Decreasing the focal ratio by a factor 
of 4 improves the sensitivity by a factor of 16. 

The beamsplitter must provide the highest possible efficiency. Uncoated pellicle beamsplitters provide rather low average 
efficiency (typically q « 0.46), although the efficiency is higher in narrow bands. Improved beamsplitters can be manufactured 
by depositing coatings on thick substrates, but substrates lacking significant absorption from 10-100 pm do not exist. Polarizing 
beamsplitters can provide ideal efficiency in a properly designed (Martin-Puplett) FTS, but such designs have yet to demonstrate 
adequate performance at wavelengths as short as 10 pm. However, bilayer pellicle beamsplitters 5 can provide nearly ideal 
performance over the desired band, as shown in Figure 2. 

The sensitivity can be increased somewhat by using larger detectors, but in practice the detector size is limited by the ability 
to design a spectrometer with sufficient throughput to fully illuminate a 10 x 10 array of detectors with a fast beam. Larger 
detectors also tend to have longer thermal time constants and thus have insufficient bandwidth for use in a scanning FTS. 

The best presently available uncooled IR detectors at high frequencies (> 100 Hz) are pyroelectric detectors. A 1 mm 2 
detector can be obtained with D* = 1.5 x 10 9 cm Hz 0 5 /W at 100 Hz, dropping to 2 x 10 8 at a kilohertz. This is still an order of 
magnitude or more less than the thermal noise limit of 1.8 x 10 10 for an uncooled detector. 3 

The increases in sensitivity that may be realized are shown in Figure 3. For both calculations we assume x = 1 .4 s, e = 1 .0, 
8 k = 0.3 cm -1 , and a scene temperature of 230 K. For the conventional technology we use the Mylar beamsplitter shown in 
Figure 2 and the 1 mm 2 pyroelectric detector described above, illuminated by an f = 2 beam. For the proposed technology 
we use the bilayer beamsplitter and a thermal noise limited detector illuminated by a parabolic feed horn. We assumed a feed 
horn with an input focal ratio of 4 and an input area of 1 mm 2 . The output focal ratio is 0.5 and illuminates a detector area of 
0.0625 mm 2 (required to conserve A£2). 





Figure 3. Comparison of predicted performance with conventional and proposed technologies described in the text. 

5. CONCLUSION 

A space-borne imaging FTS has been proposed that can provide spectra covering the wavelength range 10-100 /jm with 

0.6 cm -1 resolution, daily global coverage, a 10 km footprint, and an NEdT of better than 0.1 K from 10-50 /./in. The proposed 
instrument requires infrared detectors that are thermal noise limited at ambient temperatures. Such detectors are an order of 
magnitude more sensitive than the best detectors presently available. 
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